Entanglement of Formation of Bipartite Quantum States 
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We give an explicit tight lower bound for the entanglement of formation for arbitrary bipartite 
mixed states by using the convex hull construction of a certain function. This is achieved by 
revealing a novel connection among the entanglement of formation, the well-known Peres-Horodecki 
and realignment criteria. The bound gives a quite simple and efficiently computable way to evaluate 
quantitatively the degree of entanglement for any bipartite quantum state. 
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Quantum entangled states are used as key resources 
in quantum information processing and communication, 
such as in quantum cryptography, quantum teleporta- 
tion, dense coding, error correction and quantum com- 
putation 0]. A fundamental problem in quantum infor- 
mation theory is how to quantify the degree of entangle- 
ment in a practical and operational way |2j,|3| . One of the 
most meaningful and physically motivated measures is 
the entanglement of formation (EOF) 0, 13 , which quan- 
tifies the minimal cost needed to prepare a certain quan- 
tum state in terms of EPR pairs. Related to the EOF 
the behavior of entanglement has recently been shown 
to play important roles in quantum phase transition for 
various interacting quantum many-body systems and 
may significantly affect macroscopic properties of solids 
|{| . Moreover, it has been shown that there is a remark- 
able connection between entanglement and the capacity 
of quantum channels [|J. A quantitative evaluation of 
EOF is thus of great significance both theoretically and 
experimentally. 

Considerable efforts have been spent on deriving EOF 
or its lower bound through analytical and numerical a 
proaches, for some limited sets of mixed states 0, IT 

Among them, the most 
noteworthy results are an elegant analytical formula for 
two qubits 0, , explicit derivations for isotropic states 
0, Werner states ^l| and Gaussian states with cer- 
tain symmetries Closed-form expressions have also 
been given for special classes of high dimensional states 
jT3. flp.fl^| and rotationally symmetric states in 2(gm sys- 
tems 1 1 51 , where n is the dimension of the Hilbert space 
associated with the second subsystem. Possible lower 
bounds have been given in 0] only for 2 £g> n states. 
Notable progress has been achieved in 0, 0] in giv- 
ing analytic lower bounds [that can be optimized fur- 
ther numerically [l8j ] for the concurrence, which permits 
to furnish a lower bound of EOF for a generic mixed 
state. However, this lower bound is not explicit except 
for the case of 2 ® n system s 1 191 . For low dimensional 
systems, numerical methods [lfj can be used to estimate 
EOF. Nevertheless, they are generally time-consuming 
and often not very efficient. The notorious difficulty of 



evaluation for EOF is due to the complexity to solve a 
high dimensional optimization problem, which becomes 
a formidable task, as the dimensionality of the Hilbert 
space grows. 

In this Letter, we present the first analytical calcula- 
tion of a tight lower bound of EOF for arbitrary bipartite 
quantum states. An explicit expression for the bound is 
obtained from the convex hull of a simple function, based 
on a known result in || . This is achieved by establishing 
a key connection between EOF and two strong separa- 
bility criteria, the Peres-Horodecki criterion |2(J, |2lJ and 
the realignment criterion [2^, HH • The bound is shown to 
be exact for some special states such as isotropic states 
0, EjI and permits to provide EOF estimations for many 
bound entangled states (BES). It provides a very simple 
computable way for getting information on the actual 
value of EOF, and in particular, fills significantly the 
large gap between the nice result on the two qubits case 
and a few other existing results (mentioned above) 
for high dimensional mixed states. 

Let us first recall some useful notations. A pure m®n 
(m < n) quantum state \ip) is a normalized vector in the 
tensor product H.a®7~Lb of two Hilbert spaces Ha, Hb for 
systems A, B. The entanglement of formation is defined 
to be E(\<ip)) = S(p A ) where p A = Tr B (\ip) (ip\) is the 
reduced density matrix. Here S(pa) is the entropy 



S(p A ) = - ^/x, log 2 m = H(ji), 



(1) 



where pi are the eigenvalues of pA and p is the Schmidt 
vector (pi, p2, . . . , p m )- It is evident that E(\ip}) vanishes 
only for product states. This definition can be extended 
to mixed states p by the convex roof, 



E{p) 



(2) 



for all possible ensemble realizations p = X)i-P«l' ( /'i)('0i|) 
where pi > and J2iPi = 1- Consequently, a state p 
is separable if and only if E(p) = and hence can be 
represented as a convex combination of product states 
as p = Y^iPiPt ® Pi i where pf and pf are pure state 
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density matrices associated to the subsystems A and B, 
respectively j2j|- The measure Eq. J2J satisfies all the 
essential requirements of a good entanglement measure: 
convexity, no increase under local quantum operations 
and classical communications on average, no increase un- 
der local measurements, asymptotic continuity and other 
properties 00- 

A central idea for our approach is as follows: instead 
of the conventional method to make optimization subject 
to a large number of constraints in Eq. (J2J), we look for 
minimal admissible E{\tpi)) for a given |V>i) with some 
restrictions generated from the simple computable Peres- 
Horodecki criterion of positivity under partial transpose 
(PP T criterion) ji^, EI and the realignment criterion 
[22J, |23j . One then expects to obtain a tight lower bound 
through a convex roof construction for the pure states. 
Without loss of generality, we suppose that one has a 
general pure m® n quantum state, which can always be 
written in the standard Schmidt form 



(3) 



where ^fpi (i = 1, . . . m) are the Schmidt coefficients, 
I a;) and \bi) are the orthonormal basis in TLa and Hb, 
respectively. It can be straightforwardly verified that the 
reduced density matrices pa and ps have the same eigen- 
values p,i. It follows from Eq. Q that E(\ip)) vanishes 
only for a product state and reaches its maximum log 2 m 
for a maximally entangled state. 

Let us recollect some details on the two above men- 
tioned criteria. Peres made firstly a distinguished 
progress in the study of separability [20j by showing that 
P Ta > should be satisfied for a separable state, where 
p TA stands for a partial transpose of p with respect to 
the subsystem A. It was further shown by Horodecki et 
al. [2l| that p TA > is also sufficient for separability 
of 2 (gi 2 and 2 (g) 3 bipartite systems. In addition ||p TA || 
was shown in Refs. [5(1 Efil ] to be invariant under local 
unitary transformation (LU), where || • || stands for the 
trace norm defined by ||G|| = Tr(GG^ i ) 1 /' 2 . It is clear 
that p = Y,ijkl PikjiWifo) (ajh\ and pT^ = p ]k il in a 
suitably chosen orthonormal basis a t (i < m) and 
(k < n). Here the subscripts i and j can be regarded as 
the row and column indices for the subsystem A respec- 
tively, while k and / are such indices for the subsystem 
B. 

The realignment criterion is another powerful oper- 
ational criterion for separability given in j^, H3|. It 
demonstrates a remarkable ability to detect many BES 
[22I and even genuinely tripartite entanglement [13 • 
Recently considerable efforts have been made in propos- 
ing stronger variants and multipartite generalizations for 
this criterion [2^] . The criterion says that a realigned ver- 
sion 7Z(p) of p should satisfy ||7?.(p)|| < 1 for any separa- 
ble state p. 1Z(p) is simply defined to be TZ-(p)ij,H — Pikji 
The LU invariant property holds also for 



Die state p. , 

0SH3. 



11^0)11 E3- For the pure state of Eq. it is straight- 
forward to prove 



\p Ta \\ = 



m°)ii = (£- 



X. 



(4) 



as already shown in |22t |26| and [19| , where A varies from 
1 to m. 

Let us also review some important results in which 
we will use for the proof of our main Theorem. Tcrhal 
and Vollbrecht gave the first formula for the entangle- 
ment of formation for a class of mixed states in arbitrary 
dimension d: the isotropic states 



PF 
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where |*+) = y/T/dYli=i \ U ) and F = (^ + \Pf |* + ), 
satisfying < F < 1, is the fidelity of pF and |^ + ). They 
found that for F > 1/d, the EOF for isotropic states is 
E(pf) — co[R(F)] where R(F) is a simple function of 
F. Here the symbol "co" means the convex hull, which 
is the largest convex function that is bounded above by 
a given function. They have also presented an explicit 
expression of co[i?(F)] for d — 2,3, and conjectured its 
general form for arbitrary d. 

With the above analysis and preparation, we can for- 
mulate the main result of this Letter: 
Theorem. — For any m®n (m < n) mixed quantum state 
p, the entanglement of formation E(p) satisfies 



where 



E(p) > co[R(A)}, 

R(A) = tf 2 [7(A)] + [I - 7 (A)] log 2 (rn - f) 
7 (A) = ^[VA+ V(m-f)(m-A)] 2 , 



(6) 



(7) 



with A = max(||p TA || , ||7?.(p)||) and i?2(-) is the standard 
binary entropy function. 

Proof. — To obtain the desired lower bound, we assume 
that one has already found an optimal decomposition 
'YliiViP 1 f° r P f° achieve the infimum of E(p), where p l are 
pure state density matrices. Then E(p) = ^2 i PiE(p l ) by 
definition. For a pure state density matrix a = IV'XV'I 
with \ip) given by Eq. (0, one has ||ct Ta || = = 
(Sfeli \fp~k) 2 = A according to Eq. (gj). We would like 
firstly to find a minimal admissible H(p) for a given A. 
This minimization problem has been solved in [9j, 



7fl 

R(X) = mm{H(p) \ (jT. 



'Pk 



k=l 



A 



= H 2 [7(A)] + [1- 7(A)] log 2 (m-l), (8) 



where 



(A) - ^[VA+ V(m-I)(m-A)] 2 



(9) 
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The function R(X) here is the R(F) used in after sub- 
stitutions of d — > m, F — > A/to. It is further shown in 
that co[i2(A)] is a monotonously increasing, convex func- 
tion and satisfies co[i?(A)] < H{p) for a given A. Denote 
£ (A) = co[i?(A)], one thus has 



E( P ) = £>w) = Ef^)^E^( A< ) 

z i i 

"£(11/^11), 



> f 



(£» A< ) ^ 



f(iiwii), 



(10) 



where we have used the monotonicity and convexity prop- 
erties of £, and convexity of the trace norm ||p A \\ < 
EtPiWiP^W and \\K(p)\\ < EiPi\W)l Set A = 



max[||p 



\lZ(p)\\], we arrive at 

E{p) > £(A) = co[i?(A)] 



(11) 



which gives exactly the conclusion of Eq. JBJ . ■ 
Since co[i?(A)] is the largest convex function that is 
nowhere larger than R(X), it is optimal to give the best 
lower bound according to the relations Eqs. (JSJ), (jTUl . and 
From the general form for co[i?(A)] given in |9(, the 
following relation 

E(p) > 

0, A = l, 

H 2 [ 7 (A)\ + [l- 7 (A)]log a (m-l), Ae [l, ^% 

i25^_i)(A-TO)+l0g 2 TO, A 6 [^ 



,mj, 
(12) 

holds for to = 2,3. We have strong evidence for its cor- 
rectness for arbitrary to by verifying directly, that the 
second derivative for i?(A) with respect to A goes from 
positive to negative value with only one zero point when 
A varies from 1 to to according to the analysis of ,9| . 

The result of the above Theorem and its general ex- 
pression Eq. (|12|) provide an explicit tight lower bound 
for the EOF without the need of any numerical opti- 
mization procedure. In fact, it can be done in an entirely 
straightforward manner through the computation of the 
trace norm of a certain matrix by standard linear algebra 
packages. Some further significant features arising from 
our general result are illustrated in several examples and 
the following discussions. 

Example 1: Qubit-qudit system 

When m = 2, which corresponds to a qubit-qudit sys- 
tem, one derives easily from the Theorem and Eq. i|12|) 
that 



E{ P )> 



0. 



A=l, 



ff 2 [i(l + v /l-(A-l) 2 )], Ae [1,2]. 



(13) 



This recovers previous results obtained by the authors 
in 0], where A = max[||p TA ||, ||i?(p)||], and others in 
[29| . where A = ||p TA ||. In addition Eq. JT2J) can detect 



and give lower bounds of EOF for all entangled states of 
two qubits and of qubit-qutrit systems, since the Peres- 
Horodecki criterion is necessary and sufficient for sepa- 
rability in these cases |2l| . Furthermore, whenever there 
is a two qubit state in which A — 1 is equal to the con- 
currence defined in the bound Eq. (|T3l will give the 
exact value of EOF. For example the 2 ® 2 Werner state 
[ll| fits this condition by direct verification. The bound 
Eq. (jT3|) will be particularly useful for the study of en- 
tanglement in realistic many-body physical systems. For 
example, one usually needs to monitor entanglement dy- 
namics and distribution between a spin 1/2 particle and 
the remaining parts for a solid state system or a quantum 
computing device, and our bound can be useful in this 
context. 



Isotropic states Eq. Q were firstl y pr oposed in [2j| and 



Example 2: Isotropic states 

lypn 

further properties were analyzed in |ll| . They arise natu- 
rally in some special depolarizing channel |24j | and consti- 
tute the class of U (£>£/* invariant mixed states in d®d sys- 
tems. These states have been shown to be separable for 
F < 1/d 01 . The EOF E(p) for this class of states has 
been given in 0_by an elegant extremization procedure. 
It is derived inplE^ that \\p T F A \\ = \\TZ(p F )\\ = dF for 
F > 1/d. Thus one can directly exploit the above Theo- 
rem with A = dF to see that the bound given in Eqs. © 
and (|12fl coincides with the exact value of EOF for the 
whole class of states in Eq. (JjJJ. 

Our Theorem and the general relation of Eq. l|12f) com- 
plement a number of existing approaches to make a quite 
good estimate of entanglement for BES, benefiting from 
the powerful realignment criterion which enables one to 
detect many of the BES [HE^. 

The bound can be made even better if it is much 
easier to compute a convex-roof extended entangle- 
ment measure 

CO 

or ||7£(p)||co) than the 
EOF. The extended measures in our case are defined 
by ||p^|| co = min^j^lK^PH and \\K(p)\\ co = 
min {p . !i0l } Y.iPt\\'^(p l )\\ where p = J2iPiP\ Pi > and 
YliPi = 1- They have been studied and calculated 
for some special class of states in j^, H(J- Defining 
A = max(|jp TA || co , ||7£(p)||co), cine finds that the result 
of the Theorem is still valid, since the last inequality in 

Eq. 10 holds as £ (EiPi^) > £ ( A )- This wil1 Pro- 
vide a tighter lower bound for the EOF since generally 



and \\TZ(p)\\ co > \\1l(p)\\ follow from 



ll^llco > ||^| 

the definitions. 

Due to the nonanalytic behavior of the right hand 
side of Eq. (|12|) . it is difficult to find a specific condi- 
tion under which our bound will be exact for a general 
state. Roughly, for a20n system one necessary require- 
ment is that all the p l should be equally entangled with 
E{p l ) = R(A) in an optimal decomposition for achieving 
EOF, as in this case we demand all the inequalities in 
Eq. i|l(J[l to be changed into equalities. In higher dimen- 
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sions, it is necessary that all the p l should have equal 
EOF of R(A) when 1 < A < 4(m - l)/m while there 
should be two values of EOF, log 2 m and i?(4(m — l)/m) 
for 4(m — l)/m < A < ro, as seen from Eq. (|12[) . 

Although the work |l9j permits to furnish a lower 
bound for EOF (in fact there is no explicit formula given 
there), it is by no means optimal for general states except 
for 2 ® n systems. Instead, the procedure in 0] clearly 
imposes more restrictions than that of the present work: 
it requires firstly to give a lower bound of concurrence for 
a given A, and then obtain a possible lower bound of EOF 
from the derived concurrence bound. Thus the result in 
the present Letter is optimal as long as the parameter A 
is involved only, since we have utilized the largest convex 
function co[i?(A)] that is bounded above by R(X). 

In summary, we have determined a completely ana- 
lytic lower bound of EOF for an arbitrary bipartite mixed 
state, which characterizes optimally the quantitative be- 
havior of entanglement through the well-known Peres- 
Horodecki criterion and the realignment criterion for sep- 
arability. The procedure only involves a simple compu- 
tation of matrix eigenvalues and can be done efficiently 
with the standard linear algebra packages. Our bound 
leads to exact values of EOF for some special quantum 
states and enables one to give an easy EOF evaluation 
for many BES, a task which was extremely difficult be- 
fore. We are of the opinion that the result constitutes a 
significant bridge over the big gap between the elegant 
result of Wootters for 2 qubits, and the few existing re- 
sults mentioned before for high dimensional states. In 
this way our method can yield a powerful tool for in- 
vestigating quantitatively the character of entanglement 
for practical laboratory sources (atomic, photonic, spin 
or other carriers), and for the study of many-body sys- 
tems. Furthermore, it may provide important insights 
into realistic quantum channels and condensed matter 
systems, revealing deep connections between entangle- 
ment and macroscopic properties for the corresponding 
physical systems. 
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